Quantum magnetism in the paratacamite family: towards an ideal kagome lattice 
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We report tiSR measurements on the 5=1/2 (Cu 2+ ) paratacamite Zn :E Cu4_a : (OH)6Ci2 family. 
Despite a Weiss temperature of ~ - 300 K, the x = 1 compound is found to have no transition 
to a magnetic frozen state down to 50 mK as theoretically expected for the kagome Heisenberg 
antiferromagnet. We find that the limit between a dynamical and a partly frozen ground state 
occurs around x = 0.5. For x = 1, we discuss the relevance to a singlet picture. 

PACS numbers: 75.40.Gb, 76.75.+i 
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The idea of destabilizing a Neel state in favor of a spin 
liquid state by means of geometrical frustration, was pro- 
posed long ago by Anderson in the context of triangular 
S=l/2 antiferromagnets (AF) The underlying con- 
cept of a resonating valence bond (RVB) state, built from 
a macroscopic number of singlets, has proven since then 
to be a very rich theoretical playground and has indeed 
been advocated in the important context of hole-doped 
quantum AF such as High-T c cuprates, spin ladders, etc. 
The simpler and fundamental case of geometrically frus- 
trated AF insulators has been deeply revisited for the 
last 15 years. On a theoretical ground, there is now a 
consensus that a RVB state should be the ground state 
for 5 = 1/2 kagome @, H| or pyrochlore [J Heisenberg 
AF (HAF), but not for the more connected triangular 
lattice. The most striking results from exact diagonal- 
izations of the 5 = 1/2 kagome HAF are the fairly small 
gap value < J/20 and the existence of a continuum of 
low lying singlet excitations, between the ground state 
and the first excited triplet @, Q . 

Non-conventional dynamics were experimentally re- 
vealed in various corner sharing AF ^, 0, which 
yielded the best support to a resonating-like ground 
state but no satisfactory experimental realization has 
been yet available to check in detail the T = theo- 
retical predictions for 5=1/2 systems, e.g. the ex- 
istence of a singlet-triplet gap. Indeed, in a context 
where all perturbations such as single ion, exchange 
anisotropy or spinless defects, are relevant to induce var- 
ious marginal orderings, no material has been fulfilling 
the condition of a 5 = 1/2 perfect kagome HAF lattice. 
In [Cu 3 (titmb) 2 (OCOCH 3 ) 6 ].H 2 2 nd and 3 rd neighbor 
interactions play a dominant role [^, fiof. Volborthite is 
maybe the first 5=1/2 kagome lattice [ll| although one 
suspects spatially anisotropic nn interactions; it under- 
goes a transition to a peculiar spinglass state around 1 K 
and has persisting T — > dynamics H, HI]. As for 3D 
structures, rare earth pyrochlores are up to now recog- 



nized as perfect corner-sharing lattices but most of their 
fascinating properties are governed by the delicate bal- 
ance between magnetocrystalline anisotropy, exchange 
and, due to a large 5 value, dipolar interactions. To 
our knowledge, clinoatacamite Cu4(OH)gCl2 would be 
the first example approaching a 5 = 1/2 (Cu 2+ ) Heisen- 
berg pyrochlore lattice but it has a distorted lattice and 
possibly a ferromagnetic coupling along the c-axis. A 
recent fiSK study clearly revealed two magnetic transi- 
tions, one around 18 K and the other one around 6 K to 
a state where a weak ferromagnetic moment appears and 
where dynamics persist down to T = [lij ]. 

Quite recently, the x — 1 compound of the Zn- 
paratacamite family Zn a; Cu4_ I (OH)eCl2 has been re- 
vealed as a "structurally perfect 5 = 1/2 kagome an- 
tiferromagnet" 15j . It can be viewed as a double variant 
of the parent compound clinoatacamite (x — 0); first, 
the symmetry relaxes from monoclinic to rhombohedral 
(R3m) around x = 0.33, leading to a perfect kagome lat- 
tice in the a-b plane; then, in the c- elongated x > 0.33 
pyrochlore structure, the magnetic bridge along c-axis be- 
tween a-b kagome planes is progresively suppressed by re- 
placing the apical Cu 2+ by a non-magnetic Zn 2+ . At the 
composition x = 1, the 5=1/2 Cu 2+ ions lie at the ver- 
tices of a perfect kagome lattice. Whereas macroscopic 
susceptibility indicates that the Cu-Cu antifcrromagnetic 
interactions yield a Weiss temperature #cw ~ — 300 K, 
no signature of a transition was found down to 1.8 K. 

In this Letter, we present the first sub-Kelvin magnetic 
study of ZnCu3(OH)gCl2 using /iSR, which is very sen- 
sitive to any freezing and also perfectly suited to track 
persisting T — dynamics. We find neither evidence 
of a magnetic transition nor a drastic slowing down of 
the dynamics down to 50 mK, a very original behaviour 
in comparison to the other existing geometrically frus- 
trated magnets. We also track the crossover from a 
paramagnetic-like ground state to a T — magnetic 
frozen state from x = 1 to x = 0. 



2 



t 1 r 




time (us) 

FIG. 1: /i + polarization in x = 1 compound, either in zero 
field (H- and - D samples) and under an applied field (H- 
sample). No static internal field other than of nuclear origin 
is detected and a weak dynamical relaxation is observed. 



0.15 < x < 1.00 pure samples where prep ared using the 
hydrothermal method as described in [l5(. To make the 
deuterated analog of the x = 1 phase (D-sample) the en- 
tire reaction was carried out in a nitrogen filled glovebag, 
deuterated Copper Carbonate was prepared and 99% 
D 2 was used throughout. The purity of the samples and 
Cu to Zn ratio were verified ±1%) with powder x-ray 
diffraction and inductively coupled plasma - atomic emis- 
sion spectroscopy (ICP-AES). SQUID tests were per- 
formed and gave results in perfect agreement with those 
of ref. [15]. For x = 1, we find 9 C w = -300(20) K and 
we detect a small Curie term which would correspond to 
a maximum of 6% 5=1/2 free spins out of the total Cu 
amount. In all our x < 1 samples a ferromagnetic frozen 
component appeared which weight was found negligible 
for x = 1 as in [l^], and increased drastically for x < 0.5. 

fiSK measurements were performed on x = 0.15 - 1.0 
Zn contents, in zero field (ZF), longitudinal applied field 
(LF) or weak transverse field (wTF) - with respect to the 
muon initial polarization. Due to their positive charge, 
implanted muons come at rest in well defined sites where 
their electrostatic energy is minimized, i.e. nearby Cl _ 
and OH~. In a paramagnetic state, electronic moments 
fluctuate fast (~ 10~ 12 s) on the /iSR time scale and the 
field sensed by the muons has only a nuclear origin with 
a typical value of a few G. 

We first focus on the low-T state for the "perfect 
kagome" composition, x — 1 (T = 50 mK, fig. 1). 
Damped ZF oscillations, corresponding to a few G field 
are evident for the H-sample and were observed at all 
T, but disappear for the D-sample. They are therefore 
clearly associated with the H-nucleus. This is evidence 
for the formation of a OH-/1 complex where the muon 
binds to a hydroxyl group in the H-sample The 
signature of such a complex lies in oscillations of the po- 



larization. The characteristic frequency woh of the time 
evolution of the polarization Poh(^) can be calculated by 
modeling the dipolar interaction between the proton and 
the muon moments jl7^ . In the D-sample, the frequency 
wod is too small (nuclear moments ratio hd/^h = 0.153) 
and the oscillating pattern escapes our 0-16 /is time win- 
dow. As expected in a static case, a longitudinal field ~ 
a few 10 G, much larger than nuclear dipole fields, de- 
couples nearly completely the relaxation, except a small 
relaxing tail which is found identical in 80 G and 10 G 
data. The latter can be safely attributed to a dynami- 
cal relaxation, due to fluctuations of the electronic mo- 
ments, which only starts decoupling for much larger LF 
such as 2500 G. In addition, we do not find any change 
of the oscillations with T which evidences that there is 
no static field other than the dipolar field from nuclei. 
We can qualitatively conclude that (i) no freezing (ii) 
no drastic slowing down such as observed in kagome-like 
compounds @, H,[H, HH occur even at 50 mK. 

In order to refine this qualitative statement, the ZF 
polarization for the H-sample was fitted to the sum of 
the two expected contributions, 

(7mAoh*) 2 

P z (t) = [P OH (t)e 2 +KT cl (t)] e (- At ) (1) 

The Gaussian damping Aoh originates from nuclei sur- 
rounding the OH- \l complex. KTci(i) is the well estab- 
lished Kubo-Toyabe function, Gaussian at early times 
used here for the Cl _ site (~ 10 - 20%) where the relax- 
ation is due to surrounding Cu and H nuclei. In addi- 
tion, a small dynamical relaxation was allowed through 
the stretched exponential overall multiplying factor [18j | . 
All the static values used in the fits were fixed by a high 
statistics run at 1.5 K and found identical for higher T 
values and other x values - provided the samples are in 
a paramagnetic state. From the value of woh, we find 
i?OH = 7.6(2) G, and the field distributions are respec- 
tively A OH = 2.2(3) G and A C i = 1.7(5) G for OH" and 
Cl~ sites, typical of dipolar fields induced by nearby H, 
Cu and CI nuclei. The D-sample data can be consistently 
analyzed using -Hod = Hon X — and similar A values. 
At 50 mK we found identical Hqh values and only the 
dynamical relaxation was found to slightly increase. 

The absence of variation of the frequency of the oscilla- 
tion with T on the H-sample demonstrates that although 
6>cw ~ —300 K, the sample does not freeze into a mag- 
netic state down to the lowest T of 50 mK. Indeed, in 
the case of a magnetic freezing, the oscillation due to the 
OH-fi complex would be swamped by the relaxation due 
to electronic moments. An upper value of the electronic 
field of ~ 0.3 G corresponds to the error bars of our fits. 
From the comparison with the ^500 G field found at the 
muon site in Cu4(OH)eCl2, one obtains an upper bound 
of 6 10 _4 /is for the Cu 2+ frozen moment. 

We now turn to the discussion of the small dynami- 
cal effects observed on the long time tail of the muon- 
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FIG. 2: Left: time-evolution of the polarization for x = 1 
D-sample in ZF and LF. The slowing down observed below 
1 K is more pronounced for x — 0.66 than for x = 1 (right). 



FIG. 3: /i + polarization 19j plotted versus time for x — 
0,0.15,0.33,0.5,0.66,1.0. Lines are for fits described in the 
text. The data for x — is from 12011. 



polarization. We tracked its T-evolution by applying 
a 80 G longitudinal field, large enough to decouple the 
static nuclear fields. In the left panel of fig. 2, one can di- 
rectly notice on the D-sample raw data a small increase of 
the muon-relaxation between 1.5 and 50 mK. The values 
of A extracted from stretched exponential fits are plotted 
in the right panel of fig. 2. The modest increase of A is 
found to occur around 1 K. The stretched exponent a was 
also found to increase from 1 to 1.5 in this T range which 
reminds of the unconventional undecouplablc Gaussian 
(a=2) observed in a broad variety of geometrically frus- 
trated magnets [H, [|| . 

In the simplest model of a single time relax- 
ation process, A is related to the fluctuation fre- 
quency v and the fluctuating electronic field Hfi uct by 
A = 2 1 lH} luct v/{v 2 + 1 lH 2 LF ). From fits of the 80 G 
and 2500 G data at 50 mK (see Fig. 1), we get v ~ 
150 MHz and H } i uct ~ 18 G. This fluctuating field is 
much smaller than the static field detected in the frozen 
phase of the x = compound (~ 500 G) or the mere es- 
timate of the dipolar field created by one Cu spin at the 
oxygen site (~ 2300 G). Since muons are only sensitive to 
slowly fluctuating spins which are not paired in a singlet 
state, we can safely conclude that most of the Cu spins 
appear to be inactive at 50 mK, i.e. either in a singlet 
state, or too fast fluctuating to be sensed by the muon. 
Slowly fluctuating defects are therefore responsible for 
the small relaxation we observe. 

In order to make the bridge between the fully fluctu- 
ating ground state observed for x = 1 and the freezing 
observed for x — 0, we performed a detailed T-study of 
the evolution of the magnetic properties between these 
two compositions. In fig. 3, we plot the typical polariza- 
tions obtained at T < 1.6 K, representative of the ground 
state. The evolution is best followed starting from x = 
where the freezing of electronic moments corresponds to 
well-defined spontaneous oscillations of the polarization 



associated with the precession of the fi + in internal fields 
~ 500 G [2Cj (fig. 3, inset). For x=0.15, the oscilla- 
tions become less evident and a substantial slowly relax- 
ing component appears. The weight of the latter keeps 
increasing up to x — 0.5, at the expenses of the fast 
front end (fig. 3, main panel). The x — 0.66 sample was 
found to have a dominant dynamical character down to 
50 mK which we probed through LF measurements such 
as described for x = 1. For this composition, the slowly 
oscillating pattern of the OH-/J complex is also partially 
recovered. The slowing down of the dynamics, more pro- 
nounced than for x = 1 is clearly identified around 1 K 
(fig. 2, right panel). 

For x < 0.5, the ZF data was fitted to a superposition 
of a fast relaxing (x > 0.33) or oscillating (x < 0.15) com- 
ponent corresponding to the frozen fraction, and a para- 
magnetic term given by Eq. 1, where all the static fields 
were kept at the values found at high T consistent with 
those for x = 1 (fig. 3). The T- variation of the frozen 
fraction, determined through 20 G wTF measurements, 
is plotted in fig. 4. For x — 0.15, two transitions clearly 
occur around 18 and 6 K but, at variance with the x = 
sample, correspond to different frozen fractions. For the 
x = 0.33 and 0.5 sample, the 18 K transition disappears, 
which therefore rather seems typical of the monoclinic 
phase (x < 0.33). A very broad transition is still observed 
around the same T~6 K and with very different frozen 
fractions. The change in the shape of the transition 
curves indicates that no macroscopic phase segregation 
occurs but rather that for x < 0.5 — 0.6, islands freezing 
around 6 K coexist with paramagnetic ones which have 
slower dynamics than for x — 1. The low- T- fractions are 
represented in the inset versus x and clearly indicate a 
change of behavior from partially static to dominantly 
dynamical behavior around x=0.5-0.6. This certainly re- 
lates to the x = 0.5 percolation threshold between a limit 
of decoupled kagome planes and a case where all Cu 2+ 
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FIG. 4: T-variation of the frozen fraction for all x where a 
magnetic transition is observed. Inset: low-T fraction versus 
x. Symbols used here are the same as in fig. 3. Lines are 
guides to the eyes 

on a triangle belong to at least one fully occupied (elon- 
gated) pyrochlore pattern. Above x — 0.5 — 0.6, the 
dynamics is likely inhomogeneous due to random local 
variations in Cu occupancy of the Zn site. The extent 
of this region signals, in turn, that the physics of the 
kagome planes is only little perturbed through Cu-Cu 
interactions along the c-axis. 

ZnCu3(OH)gCl2 (x — 1) represents the first example of 
a kagome system without any 3D freezing down to tem- 
peratures as low as 1.6 10 -4 J. Its low-T behavior differs 
widely from the cases of kagome bilayers of the long stud- 
ied magnetoplumbite family [5lJ or of volborthite. We 
find no marginal spin glass state even for one order of 
magnitude smaller T / J range and the low-T relaxation 
is much smaller than the unconventional one observed in 
these spin glass-like phases [1, H, HI]. The explanation 
for this cannot be related to the absence of defects since 
this system also displays a low-T Curie tail in the SQUID 
data (~ few %/f.u.). This statement certainly needs to 
be refined by a better understanding of the nature of 
these defects in Zn- paratacamite. 

Two scenarios could be at work to explain the low- 
T dynamical properties. The first one traces back to 
the RVB initial picture or as well to any singlet ground 
state 0. If the singlet-triplet gap is large as compared 
to the T range of the data presented here, one should 
not detect any dynamics through the muon, which senses 
magnetic fields only. The remaining dynamics could be 
due to singlets which are broken by the defects and would 
yield either deconfined unpaired spins (sporadic excita- 
tions §[) as expected in a RVB model or localized un- 
paired spins, depending on the nature of the ground 
state. It is worth noting that we did not find any drastic 
change in the /iSR relaxation up to 150 K. This does not 
necessarily means the absence of a paramagnetic to sin- 
glet phase transition since in these phases the dynamics is 



either too fast for the muon or absent and the transition 
might be hard or even impossible to detect through /iSR. 
In this scenario, the modest change in the relaxation be- 
tween the x = 0.66 and 1 samples, as compared to other 
kagome compounds [1, [lj| could be well explained by an 
unchanged concentration of defects in the kagome plane 
since only the apical Zn is replaced by Cu when x is de- 
creased. The details of the interactions between defects 
and apical Cu 2+ ions might govern the low-T dynamics 
and could be responsible for the small change observed 
in the relaxation rate. The second scenario relies on the 
idea that the singlet-triplet gap is smaller than 50 mK, 
despite a 300 K J- value. The density of excitations would 
still be quite high and ZnCu3(OH) 6 Cl2 would then repro- 
duce the essential features of a " cooperative paramagnet" 
such as observed in Tb2Ti2C>7 spin liquid state Q which 
has no transition down to 70 mK. However, these sys- 
tems may only be loosely related since Tb2Ti2C»7 with 
large and Ising-like spins likely belongs to another class 
of frustrated magnets, especially if one refers to the or- 
dered spin-ice behavior found in the parent Tb2Sii2C>7 
compound [53, [23| . 

In conclusion, Zn-paratacamite is a very promising 
candidate for the S — 1/2 kagome lattice model. The 
issues of the existence of a gap and of the nature of low 
energy excitations will certainly be major avenues for fu- 
ture research in the field of geometrically frustrated mag- 
nets. The explanation of the different dynamical behav- 
iors and ground states of various kagome-like compounds 
will certainly represent a second challenge and not only 
the x = 1.0 compound is interesting but also the broad 
range of x values above 0.5-0.6 which provide a way of 
probing the physics through perturbative 3D couplings. 
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